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A measurement of elastic deeply virtual Compton scattering j*p — > 7p using e + p and 
e~p collision data recorded with the HI detector at HERA is presented. The analysed 
data sample corresponds to an integrated luminosity of 306 pb -1 , almost equally shared 
between both beam charges. The cross section is measured as a function of the virtuality 
Q 2 of the exchanged photon and the centre-of-mass energy W of the 7*p system in the 
kinematic domain 6.5 < Q 2 < 80 GeV 2 , 30 < W < 140 GeV and \t\ < 1 GeV 2 , where 
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t denotes the squared momentum transfer at the proton vertex. The cross section is deter 
mined differentially in t for different Q 2 and W values and exponential i-slope parameters 
are derived. Using e + p and e~p data samples, a beam charge asymmetry is extracted for 
the first time in the low Bjorken x kinematic domain. The observed asymmetry is attributed 
to the interference between Bethe-Heitler and deeply virtual Compton scattering processes. 
Experimental results are discussed in the context of two different models, one based on 
generalised partem distributions and one based on the dipole approach. 
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1 Introduction 



Measurements of the deep-inelastic scattering (DIS) of leptons and nucleons allow the extrac- 
tion of Parton Distribution Functions (PDFs). While these functions provide crucial input to 
perturbative Quantum Chromodynamic (QCD) calculations, they do not provide a complete 
picture of the partonic structure of nucleons. In particular, PDFs contain neither information on 
the correlations between partons nor on their transverse spatial distribution. 

Hard exclusive particle production, without excitation or dissociation of the nucleon, have 
emerged in recent years as prime candidates to address these issues [1-7]. Among them, deeply 
virtual Compton scattering (DVCS) on the proton (7*^ — ► jp) is the simplest. The DVCS 
reaction can be regarded as the elastic scattering of the virtual photon off the proton via a 
colourless exchange, producing a real photon in the final state. In the Bjorken scaling regime, 
corresponding to large virtuality Q 2 of the exchanged photon and \t\/Q 2 <C 1, where t is the 
squared momentum transfer at the proton vertex, QCD calculations assume that the exchange 
involves two partons in a colourless configuration, having different longitudinal and transverse 
momenta. These unequal momenta, or skewing, are a consequence of the mass difference 
between the incoming virtual photon and the outgoing real photon and may be interpreted in the 
context of generalised parton distributions (GPDs) or dipole amplitudes, respectively. In basic 
terms, a GPD (off-diagonal parton distribution) is the transition amplitude for removing a parton 
from the fast moving proton and reabsorbing it with a different momentum, thereby imparting 
a certain momentum transfer to the proton. In the dipole approach the virtual photon fluctuates 
into a colour singlet qq pair (or dipole) of a transverse size r ~ 1/Q, which subsequently 
undergoes hard scattering with the gluons in the proton. The t-dependence of the DVCS cross 
section carries information on the transverse momentum of partons. 

In the kinematic range of the HERA collider, where DVCS is accessed through the reaction 
e ± p — > e ± 7p [8-12], the DVCS amplitude is mainly imaginary [2], while the change of the 
amplitude with energy gives rise to a small real part. This reaction also receives a contribution 
from the purely electromagnetic Bethe-Heitler (BH) process, where the photon is emitted from 
the electron. The interference between DVCS and BH processes allows the extraction of the 
real part of the amplitude. In addition, the real part of the DVCS amplitude can be related to its 
imaginary part using dispersion relations. In the high energy limit at low momentum fraction 
x, the dispersion relations take a simple form [13] which can therefore be used for the DVCS 
process to verify the consistency between measurements of the real and imaginary parts of the 
amplitude. 

This paper presents a measurement of DVCS cross sections as a function of Q 2 and the 
Yp centre-of-mass energy W. The single differential cross section da/dt is also extracted. 
The data were recorded with the HI detector in the years 2004 to 2007, during which period 
HERA collided protons of 920 GeV energy with 27.6 GeV electrons and positrons. The total 
integrated luminosity of the data is 306 pb _1 . The data comprise 162 pb -1 recorded in e + p 
and 144 pb -1 in e~p collisions. During this HERA II running period, the electrorQ beam was 
longitudinally polarised, at a level of typically 35%. For this analysis, the periods with left- 
handed and right-handed beams are combined and the analysed data samples have a left-handed 

1 In this paper the term "electron" is used generically to refer to both electrons and positrons, unless otherwise 
stated. 
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residual polarisation of 1% and 5% for e + p and e~p collisions, respectively. Cross section 
measurements are carried out in the kinematic range 6.5 < Q 2 < 80 GeV 2 , 30 < W < 
140 GeV and \t\ < 1 GeV 2 . The range in x ~ Q 2 /W 2 of the present measurement extends 
from 5 • 10~ 4 to 10 -2 . The cross section measurements of this analysis supersede those of a 
previous HI publication [8], in which less than half of the present HERA II data was used. It is 
complementary to measurements performed at lower Q 2 using HERA I data [10]. In addition, 
using both beam charges, the beam charge asymmetry of the interference between the BH and 
DVCS processes is measured for the first time at a collider. 



2 Theoretical Framework 



In this paper, cross section measurements are compared to predictions based either on GPDs 
or on a dipole approach. At the present level of understanding, the pure GPD approach and 
dipole models, based on the proton-dipole amplitude, are not connected. However, in the low 
x domain, dipole amplitudes could be used to provide parameterisations for GPDs at a certain 
scale [14]. In this context, the DVCS process is interesting as calculations are simplified by the 
absence of an unknown vector meson wave function. The GPD model [6] used here has been 
shown to describe previous data. It is based on partial wave expansions of DVCS amplitudes and 
is a first attempt to parametrise all GPDs over the full kinematic domain. The dipole model [15], 
with a limited number of parameters, describes a large panel of low x measurements at HERA, 
from inclusive to exclusive processes. In this model, mainly using the gluon density extracted 
from fits to F 2 data, the DVCS cross section is computed using a universal dipole amplitude. 

For GPD models, a direct measurement of the real part of the DVCS amplitude is an impor- 
tant issue, as it gives an increased sensitivity to the parameterisation of the GPDs [2,6]. Indeed, 
a calculation of the real part of the DVCS amplitude requires a parametrisation of the GPDs 
over the full x range. Considering the large flexibility in the parameterisation of the GPDs, this 
is an important quantity to qualify the correct approach with GPDs. In the dipole approach, as 
the dipole amplitude refers only to the imaginary part, the magnitude of the real part can be 
predicted using a dispersion relation. 

In high energy electron-proton collisions at HERA, DVCS and BH processes contribute to 
the reaction e ± p — > e ±r yp. The BH cross section is precisely calculable in QED. Since these two 
processes have an identical final state, they interfere. The squared photon production amplitude 
is then given by 



— I-^bhI + |^4dvcs| A DVCS A BH -\- A DVCS A BH , (1) 



where A BH is the BH amplitude, A DVCS represents the DVCS amplitude and I denotes the 
interference term. In the leading twist approximation, the interference term can be written quite 
generally as a linear combination of harmonics of the azimuthal angle <fi. As defined in [2], <fi is 
the angle between the plane containing the incoming and outgoing leptons and the plane formed 
by the virtual and real photons. For an unpolarised proton beam and if only the first harmonic 
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in cos 4> and sin 0, which are dominant at low x [6], are considered, the interference term I can 
be written as 



/ oc -C [ai cos(pReA DVC s + ^Pi sin(plmA DVCS }, (2) 

where C = ±1 is the charge of the lepton beam, Pi its longitudinal polarisation and ai and a 2 
are functions of the ratio of longitudinal to transverse virtual photon flux [1-6]. Cross section 
measurements which are integrated over are not sensitive to the interference term. The mea- 
surement of the cross section asymmetry with respect to the beam charge as a function of 
allows to access the interference term. The beam charge asymmetry (BCA) of the cross section 
is defined as 

da + /d(f)-da-/d(f) 
CW) da+/d<j> + da-/d<j>' K) 

where da + /d(f) and da~/d(j) are the differential ep — > epry cross sections measured in e + p and 
e~p collisions, respectively. 

Considering the low residual polarisation of the data and the theoretical expression of a x 
and a 2 [2], a\ 3> a 2 Pi and the contribution of the sin term is neglected. Therefore, A c {4>) can 
be expressed as 

A c ((f)) =p x cos (j) = 2 A BH — — ReA 2 D ^ S . ^cos0. (4) 

\AdVCS\ + \ A BH\ 2 

The term | A dvcs 1 2 can be derived directly from the DVCS cross section measurement ctdvcs = 
\A 2 DVCS \/(lGirb), where b is the slope of the exponential t-dependence e -6 '*' of the DVCS cross 
section. As the BH amplitude is precisely known, the measured asymmetry is directly propor- 
tional to the real part of the DVCS amplitude and the ratio between real and imaginary parts 
of the DVCS amplitude, p = ReA DVC s /^A DVC s, can be extracted. This ratio p can also 
be derived using a dispersion relation [6, 16]. In the high energy limit, at low x and when the 
W dependence of the cross section is parameterised by a single term W 5< - Q \ the dispersion 
relation can be written as [13] 

p = ReA DVCS /\m.A DVCS = tan \ ^ - ) . (5) 



The ratio p can therefore be determined directly from the energy dependence of the DVCS 
cross section parameterised by S(Q 2 ). Comparison between the p values calculated from the en- 
ergy dependence of the DVCS amplitude and from its real part therefore provides an important 
consistency test of the measured BCA. 
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3 Experimental Conditions and Monte Carlo Simulation 



A detailed description of the HI detector can be found in [17]. Here, only the detector compo- 
nents relevant for the present analysis are described. HI uses a right-handed coordinate system 
with the z axis along the beam direction, the +z or "forward" direction being that of the outgo- 
ing proton beam. The polar angle 9 is defined with respect to the z axis and the pseudo-rapidity 
is given by r] = — In tan 9/2. 

The SpaCal [18], a lead scintillating fibre calorimeter, covers the backward region (153° < 
9 < 176°). Its energy resolution for electromagnetic showers is a (E)/E ~ 7.1%/y/E/GeV® 
1%. The liquid argon (LAr) calorimeter (4° < 9 < 154°) is situated inside a solenoidal 
magnet. The energy resolution for electromagnetic showers is a (E)/E ~ 11%/y/E/GeV 
as obtained from test beam measurements [19]. The main component of the central track- 
ing detector is the central jet chamber CJC (20° < 9 < 160°) which consists of two coax- 
ial cylindrical drift chambers with wires parallel to the beam direction. The measurement of 
charged particle transverse momenta is performed in the magnetic field of 1.16 T, with a res- 
olution of ap T /P T = 0.002P T /GeV © 0.015. The innermost proportional chamber CIP [20] 
(9° < 9 < 171°) is used in this analysis to complement the CJC in the backward region for 
the reconstruction of the interaction vertex. The forward muon detector (FMD) consists of a 
series of drift chambers covering the range 1.9 < r\ < 3.7. Primary particles produced at larger 
r\ can be detected indirectly in the FMD if they undergo a secondary scattering with the beam 
pipe or other adjacent material. Therefore, the FMD is used in this analysis to provide an ad- 
ditional veto against inelastic or proton dissociative events. The luminosity is determined from 
the rate of Bethe-Heitler processes measured using a calorimeter located close to the beam pipe 
at z = — 103 m in the backward direction. 

A dedicated event trigger was set up for this analysis. It is based on topological and neural 
network algorithms and uses correlations between electromagnetic energy deposits of electrons 
or photons in both the LAr and the SpaCal [21]. The combined trigger efficiency is 98%. 

Monte Carlo (MC) simulations are used to estimate the background contributions and the 
corrections for the QED radiative effects and for the finite acceptance and the resolution of the 
detectors. Elastic DVCS events in ep collisions are generated using the Monte Carlo generator 
MILOU [22], based on the cross section calculation from [23] and using a t-slope parame- 
ter b = 5.4 GeV~ 2 , as measured in this analysis (see section lOT ). The photon flux is taken 
from [24]. Inelastic DVCS events in which the proton dissociates into a baryonic system Y are 
also simulated with MILOU setting the t-slope b ine i to 1.5 GeV 2 , as determined in a dedicated 
study (see section |6\2|) . The Monte Carlo program COMPTON 2.0 [25] is used to simulate 
elastic and inelastic BH events. In the generated MC events, no interference between DVCS 
and BH processes is included. Background from diffractive meson events is simulated using the 
DIFFVM MC generator [26]. All generated events are passed through a detailed, GEANT [27] 
based simulation of the HI detector and are subject to the same reconstruction and analysis 
chain as are the data. 
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4 Event Selection 



In elastic DVCS events, the scattered electron and the photon are the only particles that are 
expected to give signals in the detector. The scattered proton escapes undetected through the 
beam pipe. The selection of the analysis event sample requires a scattered electron and a pho- 
ton identified as compact and isolated electromagnetic showers in the SpaCal and in the LAr, 
respectively. The electron candidate is required to have an energy above 15 GeV. The photon is 
required to have a transverse momentum P T above 2 GeV and a polar angle between 25° and 
145°. Events are selected if there are either no tracks at all or a single central track which is 
associated with the scattered electron. In order to reject inelastic and proton dissociation events, 
no further energy deposit in the LAr calorimeter larger than 0.8 GeV is allowed and no activity 
above the noise level should be present in the FMD. The influence of QED radiative corrections 
is reduced by the requirement that the longitudinal momentum balance E — P z be greater than 
45 GeV. Here, E denotes the energy and P z the momentum along the beam axis of all measured 
final state particles. To enhance the DVCS signal with respect to the BH contribution and to 
ensure a large acceptance, the kinematic domain is restricted to 6.5 < Q 2 < 80 GeV 2 and 
30 < W < 140 GeV. 

The reconstruction method for the kinematic variables Q 2 , x and W relies on the measured 
polar angles of the final state electron and photon (double angle method) [8]. The variable t 
is approximated by the negative square of the transverse momentum of the outgoing proton, 
computed from the vector sum of the transverse momenta of the final state photon and the 
scattered electron. The resolution of the t reconstruction varies from 0.06 at low \t\ to 0.20 GeV 2 
at high \t\. 

The selected event sample contains 2643 events in e + p and 2794 events in e~p collisions, 
respectively. Distributions of selected kinematic variables are presented in figured] for the full 
sample from e ± p collisions and compared to MC expectation normalised to the data luminosity. 
A good description of the shape and normalisation of the measured distributions is observed. 
The analysis sample contains contributions from the elastic DVCS and BH processes, as well as 
backgrounds from the BH and DVCS processes with proton dissociation, ep — > cyY, where the 
baryonic system Y of mass M Y is undetected. The sum of the latter contributes to 14 ± 4% of 
the analysis sample, as estimated from MC predictions. Backgrounds from diffractive cu and </> 
production decaying to final states with photons are estimated to be negligible in the kinematic 
range of the analysis. Contamination from processes with low multiplicity n° production was 
also investigated and found to be negligible. 



5 Cross Section and Beam Charge Asymmetry Measurements 

The full e ± p data sample is used to measure the DVCS cross section integrated over 0. The sep- 
arate e + p and e~p data samples are used to measure the beam charge asymmetry as a function 
of</>. 

The DVCS cross section, ^*p — > jp, is evaluated in each bin i at the bin centre values 
Q?, W{, U using the expression 
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(j\T°bs _ jY BH _ jyDVCS-inel-j 

O'DVCsiQi ,Wi,U) = ^DVCS-el a DVCS-el(Qi > W h U) , (6) 

i 

where N° bs is the number of data events observed in bin i. The other numbers in this equation 
are calculated using the MC simulations described in section [3j N® u denotes the number of 
BH events (elastic and inelastic) reconstructed in bin i and normalised to the data luminosity, 
^Dvcs-inei the number of i ne l as tic DVCS background events, jV° vcs ~ el the number of elastic 
DVCS events and o' 1 D y CS _ el is the theoretical 7*p — > 7p cross section used for the generation of 
DVCS events. The mean value of the acceptance, defined as the number of DVCS MC events 
reconstructed in a bin divided by the number of events generated in the same bin, is 60% over 
the whole kinematic range, for both beam charges. 

The systematic errors of the measured DVCS cross section are determined by repeating the 
analysis after applying to the MC samples appropriate variations for each error source. The main 
contribution comes from the variation of the t- slope parameter set in the elastic DVCS MC by 
±6%, as constrained by this analysis, and the 4% uncertainty of the FMD veto efficiency. These 
error sources result in an error of 10% on the measured cross section. The 20% uncertainty 
of the t-slope parameter needed to estimate the inelastic DVCS background (see section 16.21 ) 
translates into an error on the elastic cross section of 4% on average, but reaches 12% at high t. 
The modelling of BH processes by the MC simulation is controlled using the method detailed 
in [8] and is attributed an uncertainty of 3%. The uncertainties related to trigger efficiency, 
photon identification efficiency, radiative corrections and luminosity measurement are each in 
the range of 1 to 3%. The total systematic uncertainty of the cross section amounts to about 
12%. A fraction of about 85% of this error is correlated among bins. 

For the BCA measurement, the angle <f> is calculated from the reconstructed four-vectors 
of the electron and of the photon. MC studies indicate that the resolution of <p is in the range 
from 20° to 40°. The resolution of (f> is limited mainly by the resolution on the photon energy 
in the LAr and the resolution on the electron polar angle. In addition there are large migrations 
between the true and the reconstructed \<p\ from 0° to 180°, and vice versa. The asymmetry 
Ac(4>) is then determined from the differential ep — > epj cross sections da + /d<j> and da~ /d<f> 
using the formula ©. The cross sections da/dcfi are evaluated similarly to YP ~^ IP cross 
section at bin centre values 0, using the expression 



/jyobs jyBH-inel jyDVCS-ineK 

da/d^i) = 1 pVCS-cl ~ jyBH-cl) (^DVCS-el^i) + °BH-eM), ( 7 ) 



where iVf H " cl and JVf H " mel are the numbers of elastic and inelastic MC BH events, respec- 
tively, and ^Dvcs-eii^i) + ^BH-ei^i) denotes the sum of the theoretical DVCS and BH 
ep — > ep'j cross sections. Since a cos0 dependence is expected, events with < and 
4> > are combined, in order to increase the statistical significance and to remove effects 
on the asymmetry of any possible sin <p contribution from the residual lepton beam polarisation. 
The systematic error on the BCA measurement mainly arises from the part of the LAr photon 
energy scale uncertainty which is correlated between the e + p and e~p samples, estimated to be 
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0.5%. It leads to sizeable systematic errors on the measured asymmetry for close to 0° and 
180°. 

In a first step, the interference term between DVCS and BH processes, which is not known 
a priori, is not included in formula ©. In order to simulate the interference term, an asym- 
metry of the form p\ cos is added to the MC generation and passed through the full detector 
simulation and analysis chain to account for all acceptance and migration effects from true to 
reconstructed <\> values. Similarly to the data, formulae © and © are used to determine the 
reconstructed asymmetry corresponding to these MC events. To determine the value of p%, a 
X 2 minimisation is performed as a function of p± to adjust the reconstructed asymmetry in the 
MC to the measured one. MC events generated using this p 1 value are then used to correct the 
measured asymmetry for the effect of migrations. Bin by bin correction factors are determined 
from the difference between the true and the reconstructed asymmetry in the MC. 



6 Results and Interpretations 

6.1 Cross Sections and ^-dependence 

The measured DVCS cross sections as a function of W for \t\ < 1 GeV 2 and at Q 2 = 10 GeV 2 
as well as the Q 2 dependence at W = 82 GeV are displayed in figure [2] and given in table [Q 
They agree within errors with the previous measurements [8, 10-12]. The data agree also with 
models based on GPDs [6] or the dipole approach [15]. DVCS cross sections for e + p and e~p 
data are also found in good agreement with each other. As already discussed in [8], the steep 
rise of the cross section with W is an indication of the presence of a hard underlying process. 

The W dependence of the cross section for three separate bins of Q 2 is shown in figure Sa) 
and given in tabled A fit of the function W s is performed in each Q 2 bin. Figure |3jb) shows 
the obtained 5 values. It is observed that 5 is independent of Q 2 within the errors. The average 
valued 5 = 0.63 ± 0.08 ± 0.14 is in agreement with the previous measurement [8], as well as 
with the value of 5 = 0.52 ± 0.09 (stat.) measured by the ZEUS Collaboration at a lower Q 2 of 

3.2 GeV 2 [12]. 

Differential cross sections are measured as a function of t for three values of Q 2 and W 
and presented in table [3] Fits of the form da /d\t\ ~ e -6 '*', which describe the data well [8], 
are performed taking into account the statistical and correlated systematic errors. The derived 
t-slope parameters b(Q 2 ) and b(W) are displayed in figures Ufa) and (b), respectively. They 
confirm the result obtained in a previous analysis [8] and no significant variation of b with 
W is observed. Experimental results are compared with calculations from GPD and dipole 
models [6, 15]. A good agreement is obtained for both W and Q 2 dependences of the t-slopes. 
It should be noted that in the GPD model previous data of [8, 10] are used to derive the Q 2 
and W dependences of b, while no DVCS data enter in the determination of parameters of the 
dipole model. If b is parametrised as b = b + 2a' In -, with x = Q 2 /W 2 , the obtained a' value 
is compatible with and an upper limit on a' of 0.20 GeV 2 at 95% confidence level (CL) is 
derived. This value is compatible with results obtained for J/ij) exclusive electroproduction [28, 

2 Here and in all other places where results are given the first error is statistical and the second systematic. 
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29], for which the measured a' is below 0.17 GeV 2 at 95% CL. An increase of the slope 
with decreasing x (shrinkage) is therefore not observed. Such a behaviour is expected for hard 
processes and confirms that perturbative QCD can be used to describe DVCS processes. 

Using the complete analysis sample, the value of b is found to be 5.41 ± 0.14 ± 0.31 GeV -2 
at Q 2 = 10 GeV 2 . This corresponds to a total uncertainty of 6% on the (elastic) t-slope mea- 
surement for the full data sample. As in [8], this t-slope value can be converted to an average 
impact parameter of a/< r\ > = 0.64 ± 0.02 fm. It corresponds to the transverse extension of 
the parton density, dominated by sea quarks and gluons for an average value x — 1.2 ■ 10~ 3 , in 
the plane perpendicular to the direction of motion of the proton. At larger values of x (x > 0.1), 
a smaller value of a/< r\ >, dominated by the contribution of valence quarks, is estimated [4]. 

6.2 Inelastic DVCS ^-dependence 

The increased statistical precision compared to previous analyses allows a first measurement 
of the t-slope of the inelastic DVCS process. A sample of events with a signal in the FMD is 
selected. It corresponds to events with the mass of the proton dissociation system M Y in the 
range 1.4 to 10 GeV, as derived from MC studies. The contribution of inelastic DVCS events 
is extracted by subtracting the BH (elastic and inelastic) and elastic DVCS contributions, as 
estimated from the respective MC expectations. The measured differential cross section as a 
function of t is presented in figure |5] A fit of the form da/d\t\ ~ e ~ fe " ie! '*' yields b, me i = 1.53 ± 
0.26 ± 0.44 GeV -2 . In the present event sample, no indication of a dependence of bi ne i with Q 2 
or W is observed. The obtained value for bi ne i is compatible with previous determinations for 
inelastic exclusive production of p, [30] and J/ip [29]. 

6.3 Beam Charge Asymmetry 

The contributions of elastic DVCS and BH processes to the analysis sample are of similar size, 
as can be observed in figure [TJ This is a favourable situation for the beam charge asymmetry 
measurement, with a maximum sensitivity for the interference term. The measured BCA inte- 
grated over the kinematic range of the analysis and corrected for detector effects, as detailed 
in section [5] is presented in figure [6] and table @J Bins in with a size of the order of the 
experimental resolution on <p are used. 

The x 2 minimisation procedure leads to a p 1 value of p\ = 0.16 ± 0.04 ± 0.06. The result- 
ing function 0.16 cos0 is displayed in figure [6] and is seen to agree with the prediction of the 
GPD model for the first cos0 harmonic [6]. The measured asymmetry is in good agreement 
with the model prediction within experimental errors. 

As detailed in section [21 from the measured BCA and the pi value determined above, to- 
gether with the DVCS cross section, the ratio p of the real to imaginary parts of the DVCS 
amplitude can be calculated as p = 0.20 ± 0.05 ± 0.08. This is the first measurement of this 
ratio. The dispersion relation of equation © and our measurement of 8(Q 2 ) on the other hand 
leads to p = 0.25 ± 0.03 ± 0.05, in good agreement with the direct determination. While in the 
low x domain of the present measurement, the real part of the DVCS amplitude is positive, in 
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contrast, at larger x (x ~ 0.1) and lower Q 2 , a smaller and negative real part was measure 
the HERMES Collaboration [31]. 



7 Conclusion 



The elastic DVCS cross section 7*p — > 7p has been measured with the HI detector at HERA. 
The measurement is performed in the kinematic range 6.5 < Q 2 < 80 GeV 2 , 
30 < W < 140 GeV and \t\ < 1 GeV 2 . The analysis uses e + p and e'p data recorded from 
2004 to 2007, corresponding to a total integrated luminosity of 306 pb _1 , almost equally shared 
between both beam charges. The W dependence of the DVCS cross section is well described 
by a function W s . No significant variation of the exponent 5 as a function of Q 2 is observed. 
For the total sample a value 5 = 0.63 ± 0.08 ± 0.14 is determined. The steep rise of the cross 
section with W indicates a hard underlying process. The t-dependence of the cross section is 
well described by the form e~ 6 '*' with an average slope of b = 5.41 ± 0.14 ± 0.31 GeV -2 . The 
t-slopes are determined differentially in Q 2 and W and are compatible with previous observa- 
tions. The t-slope is also measured for the inelastic DVCS. The measured elastic DVCS cross 
section is compared to the predictions of two different models based on GPDs or on a dipole 
approach, respectively. Both approaches describe the data well. The use of e + p and e~p colli- 
sion data allows the measurement of the beam charge asymmetry of the interference between 
the BH and DVCS processes, for the first time at a collider. The ratio p of the real to imaginary 
part of the DVCS amplitude is then derived, directly from the measurements of the BCA and of 
the DVCS cross section to be p = 0.20 ± 0.05 ± 0.08. This ratio can also be calculated from a 
dispersion relation using only the DVCS energy dependence, leading to p = 0.25 ± 0.03 ±0.05. 
Both results are in good agreement. The GPD model considered here [6] correctly describes the 
measured BCA as well as p. The measurements presented here show that a combined analysis 
of DVCS observables, including cross section and charge asymmetry, allows the extraction of 
the real part of the DVCS amplitude and subsequently a novel understanding of the correlations 
of parton momenta in the proton. 
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Q 2 [GeV 2 ; 




ctdvcs [nb] 


8.75 


3.87 ±0.15 ±0.41 


15.5 


1.46 ±0.07 ±0.18 


25 


0.55 ±0.07 ±0.08 


55 


0.16 ±0.02 ±0.03 



W [GeV] 


ctdvcs [nb] 


45 


2.23 ±0.11 ±0.19 


70 


2.92 ±0.16 ±0.27 


90 


3.63 ±0.22 ±0.40 


110 


3.71 ±0.29 ±0.61 


130 


4.37 ±0.60 ±1.16 



Table 1: The DVCS cross section — > jp, (Jdvcs, as a function of Q 2 for W = 82 GeV and 
as a function of W for Q 2 = 10 GeV 2 , both for \t\ < 1 GeV 2 . The first errors are statistical, the 
second systematic. 





ctdvcs [nb] 


W [GeV] 


Q 2 = 8 GeV 2 


Q 2 = 15.5 GeV 2 


Q 2 = 25 GeV 2 


45 

70 

90 

110 

130 


3.06 ±0.18 ±0.25 
3.54 ± 0.29 ± 0.34 
4.93 ±0.39 ±0.52 
5.16 ±0.51 ±0.74 
5.62 ±1.34 ±1.19 


0.98 ±0.07 ±0.08 
1.46 ±0.12 ±0.12 
1.41 ±0.16 ±0.17 
1.66 ±0.23 ±0.28 
2.00 ±0.37 ±0.47 


0.31 ±0.11 ±0.05 
0.52 ±0.08 ±0.06 
0.81 ±0.13 ±0.09 
0.63 ±0.17 ±0.15 
0.80 ±0.26 ±0.29 


5 


0.61 ±0.10 ±0.15 


0.61 ±0.13 ±0.15 


0.90 ±0.36 ±0.27 



Table 2: The DVCS cross section 7*p — > 7p, a DVC s, as a function of W for three Q 2 values 
and for \t\ < 1 GeV 2 . The values of S(Q 2 ) obtained from fits of the form W 5 are given. The 
first errors are statistical, the second systematic. 
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da DV cs/d\t\ [nb/GeV 2 ] 


W = 82 GeV 




t 




GeV 2 ] 


Q2 = g GeV 2 


Q 2 = 15.5 GeV 2 


Q 2 = 25 GeV 2 


0.10 
0.30 
0.50 
0.80 


13.3 ±0.80 ±1.73 
4.82 ±0.32 ±0.50 
1.26 ±0.14 ±0.18 
0.21 ±0.03 ±0.04 


4.33 ± 0.35 ± 0.65 
1.24 ±0.13 ±0.16 
0.45 ±0.06 ±0.05 
0.10 ±0.01 ±0.02 


1.68 ±0.31 ±0.42 
0.49 ±0.10 ±0.08 
0.18 ±0.04 ±0.03 
0.05 ±0.01 ±0.01 


b [GeV- 2 ] 


5.87 ±0.20 ±0.32 


5.45 ±0.20 ±0.29 


5.10 ±0.38 ±0.37 








Q 2 = 10 GeV 2 




t 




GeV 2 ] 


W = 40 GeV 


W = 70 GeV 


W = 100 GeV 


0.10 
0.30 
0.50 
0.80 


4.77 ± 0.50 ± 0.49 
1.62 ±0.23 ±0.18 
0.69 ±0.11 ±0.07 
0.10 ±0.02 ±0.01 


7.81 ±0.51 ±0.85 
2.88 ±0.22 ±0.28 
0.91 ±0.10 ±0.10 
0.16 ±0.02 ±0.02 


11.0 ±0.85 ±2.23 
3.71 ±0.31 ±0.49 
1.18 ±0.13 ±0.16 
0.24 ±0.03 ±0.04 


b [GeV~ 2 ] 


5.38 ±0.30 ±0.23 


5.49 ±0.19 ±0.26 


5.49 ±0.20 ±0.35 



Table 3: The DVCS cross section j*p — > jp, differential in t, da DVC s/dt, for three values of 
Q 2 at W = 82 GeV, and for three values of W at Q 2 = 10 GeV 2 . Results for the corresponding 
t-slope parameters b are given. The first errors are statistical, the second systematic. 



[deg.] 


Ac(<P) 


10 


0.326 


± 0.086 


±0.180 


35 


0.119 


± 0.076 


± 0.090 


70 


-0.039 


± 0.080 


± 0.030 


110 


0.035 


± 0.092 


± 0.028 


145 


-0.234 


± 0.079 


± 0.076 


170 


-0.210 


± 0.075 


±0.169 



Table 4: The DVCS beam charge asymmetry A c (4>) as a function of and integrated over the 
kinematic range 6.5 < Q 2 < 80 GeV 2 , 30 < W < 140 GeV and \t\ < 1 GeV 2 . The first errors 
are statistical, the second systematic. 
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Figure 1 : Distributions of the energy (a) and polar angle (b) of the scattered electron, the energy 
(c) and polar angle (d) of the photon, the </> azimuthal angle between the plane of incoming and 
outgoing lepton and the plane of virtual and real photon [2] (e) and the proton four momentum 
transfer squared \t\ (f). The data correspond to the full e ± p sample and are compared to Monte 
Carlo expectations for elastic DVCS, elastic and inelastic BH and inelastic DVCS. All Monte 
Carlo simulations are normalised according to the luminosity of the data. The open histogram 
shows the total prediction and the shaded band its estimated uncertainty. 
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Figure 2: The DVCS cross section YP "~ > IP as a function of Q 2 at W = 82 GeV (a) and as 
a function of W at Q 2 = 10 GeV 2 (b). The results from the previous HI [10] and ZEUS [12] 
publications based on HERA I data are also displayed. ZEUS measurements are propagated 
from W = 104 GeV to 82 GeV using a W dependence W ^ 2 . The inner error bars represent 
the statistical errors, the outer error bars the statistical and systematic errors added in quadrature. 
The dashed line represents the prediction of the GPD model [6] and the solid line the prediction 
of the dipole model [15]. 
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Figure 3: The DVCS cross section YP ~> IP as a function of W at three values of Q 2 (a). The 
solid lines represent the results of fits of the form W 5 . The fitted values of 5(Q 2 ) are shown in 
(b) together with the values obtained using HERA I data [10]. The inner error bars represent the 
statistical errors, the outer error bars the statistical and systematic errors added in quadrature. 
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Figure 4: The fitted t-slope parameters b(Q 2 ) are shown in (a) together with the t-slope pa- 
rameters from the previous HI [10] and ZEUS [12] publications based on HERA I data. In 
(b) the fitted t-slope parameters b(W) are shown. The inner error bars represent the statistical 
errors and the outer error bars the statistical and systematic errors added in quadrature. The 
dashed line represents the prediction of the GPD model [6] and the solid line the prediction of 
the dipole model [15]. 
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Figure 5: The inelastic DVCS cross section differential in t at W = 82 GeV and Q 2 = 10 GeV 2 
and for events with 1.4 < My < 10 GeV The inner error bars represent the statistical errors, 
the outer error bars the statistical and systematic errors added in quadrature. 
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Figure 6: Beam charge asymmetry as a function of the angle <\> as defined in [2], integrated over 
the kinematic range of the analysis. The inner error bars represent the statistical errors, the outer 
error bars the statistical and systematic errors added in quadrature. The function 0.16 cos0 is 
also shown (solid line), together with the GPD model prediction (dashed line). 
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